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We consider adiabatic pumping of electrons through a quantum dot. There are two ways to operate 
the pump: to create a dc current / or to create a dc voltage V . We demonstrate that, for very slow 
pumping, I and V are not simply related via the dc conductance G as / = VG. For the case of a 
chaotic quantum dot, we consider the statistical distribution of VG — I. Results are presented for 
the limiting cases of a dot with single channel and with multichannel point contacts. 

PACS numbers: 73.23.-b, 72.10.Bg 



I. INTRODUCTION 

An electron pump is a device that converts a periodic 
variation of its characteristics into a time-independent 
electric current.El Such characteristics can be "macro- 
scopic" , like the charge on the device or the conductance 
of point contacts, or "microscopic" , such as the location 
of a scatterer or the magnetic flux threading the sample. 
When there are two characteristics of the device that can 
be varied harmonically with a frequency u and phase dif- 
ference <j), pumping of electrons already occurs in the adi- 
abatic limit ui — > 0. In that case the pumped current is 
proportional to u> sin (f>, and changes sign when the phase 
relationship between the parameters is reversed. Adia- 
batic electron pumps have been realized experimentally 
in (arrays of) Coulomb blockaded quantum dots, using 
the voltages on plunger gates, and/or theJ-ransparencies 
of the contacts as pumping parameters.Ba 

In this paper, we consider an adiabatic electron pump 
that consist of a semiconductor quantum dot coupled 
to two electron reservoirs by means of ballistic point 
contacts.B Variation of two gate voltages allow for small 
changes of the shape of the dot, and thus for the flow, 
of a dc current. Following a proposal by Spivak et al.p 
such a device has been built and investigated by Switkes 
et aln The device of Ref. [?] is referred to as an "adia- 
batic quantum pump" , because the variation of the gate 
voltages predominantly affects the quantum interference 
of the electrons in the quantum dot, not their classical 
trajectories. An important property of the quantum dot 
used in the experiment of Ref. |^ is that its classical dy- 
namics is chaotic. As a result, the magnitude and the sign 
of the expected pumped current I are subject to meso- 
scopic fluctuations. Since these fluctuations are large, the 
mean (7) and variance {I 2 ) are insufficient to describe the 
ensemble, and one needs to know the entire probability 
distribution P(I). 

Theoretical analysis has .focused on the dc current I 
pumped through the dotjBllj 

I = ^ J dtl(t), (1) 



T = 2tt/u> being the period of the pumping cycle. How- 
ever, experimentally, the preferred measurement is that 
of the dc voltage V that the electron pump generates]!! 

V = ^fdtV(t). (2) 

Naively, one might expect that, for small pumping am- 
plitudes, I and V are related via the dot's conductance G 
as I = VG. However, as we show in this paper, this naive 
"Ohm's law" is not always true, depending on whether 
the pumping frequency u is small or large compared to 
the charge-relaxation rate 7 of the reservoirs. (For adia- 
batic pumping, u> must always be small compared to the 
charge relaxation rate of the quantum dot.) For u> <C 7 
and if the number N of propagating channels in the point 
contacts between the quantum dot and the electron reser- 
voirs is small, the pumped current I and the difference 
D = VG — I can actually be of comparable magnitude. 

A qualitative explanation why the pumped voltage and 
current are not related via the simple relation VG — I in 
the limit u> <C 7 follows from the observation that both 
the pumped current and the pumped voltage have dc and 
ac components. For slow pumping, the electron pump 
generates a bias voltage that counteracts both the dc and 
ac currents generated in the dot. Since the conductance 
G itself also varies in time, V ac G has a dc component. It 
is this additional rectified dc component of the current 
that is responsible for the difference between VG and I 
for an adiabatic electron pump. When pumping is faster 
than the charge relaxation rate of the reservoirs, no ac 
bias voltage is generated to balance the ac current, and 
the difference between VG and I disappears. 

The purpose of this paper is to find the distribu- 
tion of the difference VG — I between pumped current 
and pumped voltage for adiabatic pumping of electrons 
through a chaotic quantum dot and to compare it to the 
distributions of I and V. For a chaotic dot, these dis- 
tributions have a universal form, independent of details 
of the pumping mechanism or the shape of the quan- 
tum dot. In section |H| we use the scattering approach to 
present a quantitative theory for the pumped voltage V, 
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FIG. 1. Schematic of a chaotic quantum dot whose shape 
can be changed by varying gate voltages X\ and Xi (left). In 
one cycle Xi and X2 trace out a contour in (Xi,Xs) space 
(right). 
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FIG. 2. Equivalent circuit for a measurement of the 
pumped voltage through the dot (left) and of the pumped 
current (right). Ci and C are geometrical capacitances of the 
i— th reservoir and the dot. 



and the difference D — VG — /. 
evaluate the distribution of VG 



In section III we then 
- / for an ensemble of 



chaotic quantum dots, using random matrix theory. We 
conclude in Sec. IV. 



II. PUMPED CURRENT AND VOLTAGE 

We consider a quantum dot coupled to two electron 
reservoirs via point contacts, see Fig. [I]. The shape of the 
quantum dot is varied periodically by variation of two 
gate voltages, represented by dimensionless parameters 
Xi and X 2 . Alternatively, X\ or X 2 can represent the 
value of an applied magnetic field, or any other parameter 
that characterizes the quantum dot. 

As discussed in the introduction, the electron pump 
can be characterized experimentally via a direct mea- 
surement of the pumped current, or via a measurement 
of the voltage V(t) between the two reservoirs generated 
as a result of the pumping of charge through the quan- 
tum dot. A formula for the dc current has been derived in 
Refs. ||,[| To derive a formula for the dc voltage V, we in- 
troduce a simple model for the quantum dot and the two 
electron reservoirs, see Fig. g. The dot and the reservoirs 
1, 2 are connected to a screeni ng gat e via capacitances 
C and Ci, C 2 - Following Refs. [m|Ji^ , we introduce the 
emissivity e8q(m)/6Xj, which is the charge that exits the 
dot through point contact m (m = 1,2) when the param- 
eter Xj (j = 1,2) is changed adiabatically by an amount 
SXj. Then the total current flowing through contacts 1 
and 2 reads 



Ji(i)=e£ 



2 



Sg(l) dX, 
SX, dt 



[V 1 (t)-V 2 (t))G, 



^)=e£« + [F 2 (*)-r l( MK, 



(3a) 
(3b) 



where G is the dc cconductance of the quantum dot. Here 
we assume that all variations are made slowly on the scale 
of the dwell time of the quantum dot. 



When current is measured, the voltages of the two 
reservoirs are equal, Vi(t) = V 2 (t). Hence the dc current 
/ can be found from integration of Eq. (|3a|) or (^) .H Us- 
ing Stokes' theorem, / can be rewritten as an integral 
over the surface area S enclosed by the contour of the 
parameters Xi and X 2 in the (Xi, X 2 ) plane, 



where 



1= — I dX 1 dX 2 i(X 1 ,X 2 ), 



d dq(l) d dq{l) 



(4a) 



dX 2 dXi dX 1 dX 2 

d dq{2) d dq{2) 
~ dX 2 dXi + dX 1 dX 2 ' 



(4b) 



When voltage is measured, the result depends on whether 
variation of the parameters Xi and X 2 is fast or slow 
compared to the charge relaxation rates 71,2 ~ GjC\. 2 of 
the reservoirs. If the variation is fast compared to 7^2 
(but still slow compared to the charge relaxation rate 
of the quantum dot), the voltage difference V\ — V 2 is 
essentially time independent and takes the value 

V = T/G, (5) 

where G is the conductance averaged over one cycle, 

G = ^£ dtG(X 1 (t),X 2 (t)). (6) 

In the opposite limit u> <^ 71,2, which is the case 
we'll consider in the remainder of the paper, one has 
I\{t) = T]I 2 (t), where 77 = C\/C 2 is a numerical coeffi- 
cient describing the capacitive division between the two 
reservoirs. Combining this with Eq. (^), we find 



Vi(t)-V 2 (t) = 



Gii + ^^ysx 



dX, 
dt 



5q(2) dXi 

1 sXi ~dT 



(7) 



The dc voltage V is then found by integration over the 
period T, with the result 



2 



T f> 

V = ^- e dX 1 dX 2 v(X 1 ,X 2 ), 



d Adq(l) 



d 



dX 2 y g dX 1 ' dX^g dX : 
d f ldq(2). d ,1 dq(2) 



I + 77 



dX 2 g dXi ' dX 1 y g dX 2 



(8b) 



where g = hG/2e 2 is the dimensionless conductance. For 
small harmonic variations of the parameters, Xi(t) — 
6X1 cos(wt) and X 2 (t) = 5X 2 cos(u}t+4>), the integrations 
are trivial, and one finds 



/ = — eui5Xi5X 2 i sin < 

for the pumped current, and 
hid 

V = — 5X x 5X 2 vsiTid 
4e 



(9) 



(10) 



for the pumped voltage. 

If the conductance G were constant in a cycle, i.e., G 
would not depend on Xi and X 2 , Eqs. ([It]) and (Sb)-jlC|) 
give the identity VG = I: the pumped current measured 
at zero bias and the pumped voltage measured at zero 
current are related by Ohm's law. However, in general G 
does depend on X\ and X 2 , and this "Ohm's law" does 
not hold. The deviation is described by the difference 

D = VG-I. (11) 

For small 8X X and 6X 2 , we find from Eqs. @), Q, 

(12a) 



1 



D = — eu}6Xi6X 2 dsm < 



d = 



V 



1 dg dq(l) 1 dg dq{l) 
gdX x dX 2 gdX 2 dX x 
1 3g dq(2) 1 8g dq(2) 
gdXi dX 2 gdX 2 dXi 



(12b) 



The derivatives to X\ and X 2 appearing in the above 
formulae are derivatives taken at constant values of the 
electro-chemical potential /x of the reservoirs. These are 
not necessarily equal to derivatives taken at a constant 
value of the (self-consistenttj) electrostatic potential V sc 
inside the quantum dot. Changes of V^c occur in a pump- 
ing cycle, because the total charge on the dot may vary 
during the pumping cycle. For technical reasons, it is 
preferred to treat V^ c , or, equivalently, the kinetic en- 
ergy E = fi — V sc , as an independent parameter, and 
to take derivatives at constant E. The above equations 
for pumped voltage and current can be rewritten using 
derivatives at constant E if we substitute the paramet- 
ric derivatives d/dX l and d/dX 2 in Eqs. ©, @), and 



(121) byt 



d d dE d 

~dX ^ * ~dX ^ dXdE' 
dE _ dg/dXj 
OX, ~ ~ C/2e 2 + dq/dE' 



(13a) 
(13b) 



Here C is the capacitance of the quantum dot. In Eq. 
(13a), we abbreviated 



dg _ flg(l) dq(2) 
dXj dXj dXj ' 



(14) 



For a realistic quantum dot, the charging energy e 2 /C is 
much larger than the mean level spacing A. In that limit, 
one finds that the dot charge remains constant during 
the pumping cycle, I\{t) = —I 2 {t) for all time. As a 
consequence, the pumped voltage V and the difference 
D = VG — I lose their dependence on the capacitive 
division 77. 

In the absence of inelastic processes and for low tem- 
peratures (temperature T below the mean level spacing 
A in the quantum dot), the emissivities dq{m) / dXj and 
the conductance G can be expressed in terms of the scat- 
tering matrix S of the quantum dot and its derivatives 
to X\, X 2 , and E. The matrix S has dimension 2N, 
where N is the number of propagating channels in each 
point contact; it is unitary (unitary symmetric) in the 
presence (absence) of a time-reversal symmetry breaking 
magnetic field. The derivatives of S are_parameterized 
via hermitian matrices R, Rj, defined asEj 



2tt dE 



dX, 



Then the emissivities dq(m)/dE, dq(m)/dX, m = 1,2, 
are given byliJ 



dq(m) 
dE 

dq(m) 
dXi 



^■Re tr P m R, 
1 

—Re tr P m Rj. 
Zn 



(15a) 
(15b) 



Here Pi = 1 — P 2 is a diagonal matrix with elements 
(-^1 )n = 1 if 3 — 171 ana - zer0 otherwise. The dc conduc- 
tance g is given by the Landauer formula, 
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= tr5 t Pi5P 2 - 



(16) 



In the next section, we shall study the distribution of 
the dimensionless difference d = vg — i for the case of a 
chaotic quantum dot. 



III. DISTRIBUTION OF D = VG — I 

For an ensemble of chaotic quantum dots, the statis- 
tical distribution of the scatteringp«atrix and its deriva- 
tives is known from the literature £3 It takes its simplest 
form when the derative of S is not parameterized by the 
matrices Rj, R, but by the symmetrized derivatives Q, 
Qj, 

Q = S- 1/2 RS 1/2 , Q 3 = S- 1/2 RjS 1/2 . 

In the prensence of time-reversal symmetry (labeled by 
the Dyson parameter = 1), the matrices Q, Q\, and 



3 



Qi are real symmetric. When time-reversal symmetry is 
broken by a magnetic field ((3 = 2), Q, and Qi, and Q 2 
are hermitian. For ideal contacts, the joint distribution 
P(S,Q,Qi,Q 2 ) readsH 



Pa 1 



2e^ 
CA' 



tr Q (det Q) 



-N/2-2(0N+2-0) 



x exp 
x exp 



e(Q) 



(17) 



^tr((Q- 1 Q 1 ) 2 + (Q- 1 Q 2 ) 2 



where 6(Q) = 1 if all eigenvalues of Q are positive and 
Q(Q) = otherwise. The capacitance C appears in the 
distribution because the ensemble is obtained by sweep- 
ing an external gate voltage, not the self-consistent en- 
ergy e)b 

To find the distribution P(d), we first integrate over 
Qi and Q 2 at fixed S and Q, and then over S and Q. 
The first integration can be done analytically, since, for 
fixed Q, Qi and Q 2 are Gaussian random matrices, see 
Eq. (|l7j). The result of this integration takes a simple 
form, 



P(d) 



J_ p -\d\/<r 

2a 



(18) 



S,Q 



where the brackets indicate the average over S and Q 
that remains to be done. In equation ([l8]) a is positive 
function of S and Q, given by 



[tvA 2 



0/3,1 



tr (PRSPR T S^- 



PRPR)] trP 2 - (trAB)' 



(19) 



where we abbreviated 



A = R[P- 



trPP 



B = R\ A- 



CA/2e 2 + trP 

tvKR 
CA/2e 2 + tri? 



and P = (P 1 -7 7 P 2 )/(l + r ? ), A = 2 (P 1 5P 2 5t-5P 2 5tp 1 ). 

For the remaining integrations over S and Q we con- 
sider two limiting cases: multichannel point contacts 
(N 1) and single channel point contacts (N = 1). 



A. Single-channel contacts 

For N = 1 the remaining number of variables is small, 
and can be integrated over numerically, using the distri- 
bution (0). Results for the (physically relevant) limit 
CA -c e 1 ^ are shown in Fig. ||. In this limit, the distri- 
bution P(d) does not depend on the capacitive division 
7] between the reservoirs. The results for larger values of 




FIG. 3. Distribution of the normalized difference 

d = vg — i between pumped voltage v and pumped current i 
for single-channel point contacts, for the physically relevant 
limit CA <C e 2 (main figure). The opposite limit CA 3> e 2 is 
shown in the inset for the case 77 = 0. Presence (absence) 
of time-reversal symmetry /3 = 1 (/3 = 2) is shown solid 
(dashed). 



C are not very different from those shown in Fig. [| This 
is illustrated in the inset of Fig. where we have shown 
the distribution P(d) for the case of large capacitance 
CA 3> e 2 and asymmetric reservoirs (capacitive division 
7] = 0). For d close to zero, the distribution shows a cusp 
with a logarithmically divergent derivative at d = 0. For 
d ^> 1, the distribution P(d) has power-law tails. For 
CA < e 2 these areB 



P(d) 



J" 3 , (3 = 1, 

tr 3 log d, p = 2. 



(20) 



The tails of the distribution correspond to samples with 
an anomalously large eigenvalue of £X, corresponding to 
an anomalously large dwell time td Ji3 a value of d in the 
tail of the contribution typically corresponds to a dwell 
time td ~ Tnd 2 ^ , th = h/A being the Heisenberg time. 
Since configurations with anomalously large dwell times 
are more sensitive to dephasing or thermal smearing, such 
perturbations will truncate the tails for d > {t^/th)^^ 2 
or d > {t h T)-' 3 / 2 . 

For comparison, we have also calculated the distribu- 
tion of the pumped current i and voltage v. For single- 
channel contacts, the distribution again takes the form 
( |l8| ) , with a replaced by <7j and a v , respectively. Expres- 
sions for <7i and a v can be found in the appendix. The 
resulting distributions are shown in Fig. §. (The dist ri- 
bution of the current was calculated previously in Ref.Jj.) 
The main conclusion upon comparison of Figs. [| and |]is 
that, for single-channel point contacts, the distributions 
of i and of d = vg — i have comparable widths. Hence, for 
N = 1, deviations from "Ohm's law" , as characterized by 
d — vg — i are of the same order as the pumped current 
i itself. 
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FIG. 4. Distribution of normalized pumped current i and 
(inset) voltage v for single-channel contacts. Presence (ab- 
sence) of time-reversal symmetry /3 = 1 (/3 = 2) is shown 
solid (dashed). 



B. Multi-channel contacts 

The distribution of d in the limit N ^> 1 can be di- 
rectly obtained from Eq. (|l8|). The integration over the 
unitary matrix U of eigenvectors of R and over the scat- 
tering matrix S is performed using the method of Ref. 
p0| . For the remaining integration over the eigenvalues 
Ti, i = 1, . . . , 2N, of the matrix R it is sufficient to know 
their density,E3 



2N 



P(r) = - r)> = —V(t + -t)(t-t^), (21) 

i — < 7TT" 
3=1 



where r± = (3 ± \/8) /27V. This gives the result 

2 / a \ 2 



(a 2 ) = 



/3N 4 



0\l + r)J \A + 2e 2 /C 



A 



(22) 



Calculation of higher moments of a shows that fluctua- 
tions of a are small compared to the average as N — > oo. 
Hence, we can conclude that the distribution of d in 
multi-channel limit is of the Poissonian form (|l8|), with 
a 2 given by Eq. @. 

This is to be contrasted with the distribution of the 
pumped current i, which is Gaussian,fpi; large N, with 
zero mean and with root mean squareErli3 



nN 2 ' 



(24) 



IV. CONCLUSION 

In summary, for adiabatic pumping of electrons 
through a chaotic quantum dot, we have derived ex- 
pressions of the pumped current I (in case of a cur- 
rent measurement) or the pumped voltage V (in case of 
a voltage measurement) in terms of the scattering ma- 
trix of the quantum dot. Pumped current and voltage 
are not simply related by the dot's conductance G. We 
have calculated the distribution of the difference VG — I 
for small pumping amplitudes, which is universal for a 
chaotic quantum dot. If the number N of propagating 
channels in the contacts between the quantum dot and 
the reservoirs is one, I and VG — I can be of the same 
size; if N ^> 1, VG — I is typically a factor N smaller 
than /. Our results are valid in the limit of slow pump- 
ing, where the pumping frequency u is much smaller than 
the charge relaxation rate 7 of the reservoirs. If u> 3> 7, 
the difference VG — I is suppressed. 

The results obtained here are important in view of the 
interpretation of the experiment of Ref. |7|. The obser- 
vations of that experiment can also be explained if the 
observed dc voltage is the result of rectification of ac dis- 
placement currents generated by the time-dependent gate 
voltages that should drive the electron pump.H There- 
fore, it is important to identify signatures that distin- 
guish adiabatic pumping from mere rectification of dis- 
placement currents. For the case of a current measure- 
ment, two such signatures are the magnetic field ¥yfp-i 
metry and the typical size of the pumped current Jl^'EHl 
Except for the case N — 1 , our results allow to translate 
these signatures to a voltage measurement as well.cj Fur- 
ther, the relation between pumped voltage and pumped 
current provides a third signature of a adiabatic pump- 
ing: For few-channel point contacts, /, VG, and the dif- 
ference VG — I are all random and of comparable mag- 
nitude for a quantum pump, while, if the dc signal is due 
to rectification, there is a fixed relationship / o^-G 2 V, 
the proportionality constant being non-universalE3 

We thank C. Marcus for discussions. This work was 
supported by the NSF under grant no. DMR-0086509 
and by the Sloan foundation. 



72x1/2 = 



1 

ttTV' 



(23) 



Hence we conclude that, for large N, typically J is a fac- 
tor ~ N smaller than the pumped current, and can be 
neglected in a measurement. Hence, in the limit N — > 00 
the expectation of "Ohm's law" vg — i holds, and one 
readily concludes that the pumped voltage v has a Gaus- 
sian distribution with zero mean and with root mean 
square 



APPENDIX A: 

Below we present results for the "intermediate" distri- 
butions of the normalized current i, voltage v, and of the 
difference d = vg — i after integration over the matrices 
Qi and Q2 at fixed S and Q, but before integration over 
S and Q, for the case of a chaotic quantum dot with two 
single-channel point contacts. For that case, the distri- 
bution of d is given by 
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P(d) 



I— 

\2<7 d 



D -\d\/a d 



(Al) 



where the brackets (. . -)s,q denote the remaining aver- 
age over the scattering matrix S and the symmetrized 
time-delay matrix Q. The distributions for i, v have the 
same form with ad replaced by <jj and a v , respectively. 
However, we should note that, unlike for the difference 
d, the form ( |A1[ ) does not hold for the distributions of i 
and v when N > 1. 

Below we list (statistical) expressions for ad, ai and 
a v for N = 1 for the case CA <C e 2 . We introduce 
the eigenvalues T\ , r 2 of the normalized time-delay matrix 
R. Their distribution can be found in Ref. Further, 
we introduce two independent random variables t (uni- 
formly distributed between and 1) and <j) (uniformly 
distributed between and 2-7r) that arise from the ran- 
domly distributed eigenvectors of R and the phases of 
the scattering matrix S. Finally, the equations for ad, at, 
and a v contain the dimcnsionless conductance g € [0, 1], 
which has distribution P(g) = {f3/2)g- 1+p / 2 . We then 
find 

2 « (nr 2 ) 3 (l- 5 ) 2 
a d,/3=l = 25 °7 



(n + t 2 ) 2 g 



af*-i = 256 
■2 



(rir 2 ) 3 



9, 



(rir 2 ) 3 1 



(ti+t 2 )V' 



1 - g ( 8tit 2 



g \n + t 2 



x (nr 2 + t(l - t)(n - t 2 ) 2 sin 2 



4/3=2 = (4rir 2 ) 2 1 - 4i(l - t) 



Tl ~ T 2 
Tl + T 2 , 



2 1 / 8TiT 2 . r 2 

^,/J=a = ^3 {^T^ 2 I L T1T2 + ( Tl - T2 ) 
x 



(2^/(1 - gW - *) sin ^ + (1 - 2t)Vff) ' 
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since S{B) = S T {-B) under reversal of a magnetic field 
B, this equation results in v(—B) = v(B). There is no such 
magnetic field symmetry of the pumped voltage for N > 1. 
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